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The consequences of localized, classical magnetic moments in superconductors are explored and their 
effect on the spectral properties of the intragap bound states is studied. Above a critical moment, 
a localized quasiparticle excitation in an s-wave superconductor is spontaneously created near a 
magnetic impurity, inducing a zero-temperature quantum transition. In this transition, the spin 
quantum number of the ground state changes from zero to i, while the total charge remains the 
same. In contrast, the spin-unpolarized ground state of a d-wave superconductor is found to be stable 
for any value of the magnetic moment when the normal-state energy spectrum possesses particle-hole 
symmetry. The effect of impurity scattering on the quasiparticle states is interpreted in the spirit of 
relevant symmetries of the clean superconductor. The results obtained by the non-self-consistent (T 
matrix) and the self-consistent mean-field approximations are compared and qualitative agreement 
between the two schemes is found in the regime where the coherence length is longer than the Fermi 
length. 



PACS numbers: 74. 25. Ib, 71.27.+a, 61.16.Ch 

I. INTRODUCTION 

One of the most interesting phenomena in supercon- 
ductors is their response to defects. They can be applied 
as experimental probes in studying the properties of the 
superconducting state. The simplest example of pair- 
breaking defects is magnetic impurities in singlet super- 
conductors, where magnetic scattering disrupts pairing 
in the singlet channel. Pair-breaking defects are known 
to lead to the formation of bound quasiparticle states in 
conventional nodeless superconductors [1-3] and virtual 
bound states in gapless superconductors [4-7]. These 
states are localized in the neighborhood of the defects 
and can have a very anisotropic structure, depending on 
the energy-gap structure [4-7]. Quasiparticle states, lo- 
calized near impurities, carry spin and are expected to 
modify the ground-state properties of the superconduc- 
tors as well. Because experimental techniques are be- 
coming increasingly capable of analyzing the effect of 
impurities on the superconductivity, it is also interest- 
ing to explore theoretically the implications of impurity 
scattering together with the symmetries present in the 
problem. Indeed, the progression from tunnel-junction 
spectroscopy [8] to scanning-tunneling microscopy [9-11] 
provides a strong evidence of inhomogeneous quasiparti- 
cle states induced by defects in superconductors. 

In this paper, we revisit the problem of a localized mag- 
netic moment interacting with a superconductor. A re- 
markable aspect of this interaction is the first-order zero- 
temperature transition which takes place in an s-wave 
superconductor as a function of the "magnetic moment" , 



w = JS/2, where S is the local impurity spin and J is the 
exchange coupling. In this transition, the spin quantum 
number s of the electronic ground state of the supercon- 
ductor l^o}™ changes from zero for a subcritical moment 
w < w c to i for w > w c . The total spin becomes S ± ^ 
depending on the sign of the exchange interaction J. The 
first to point out the phase transition was Sakurai [12] 
who showed that the transition corresponds to a level 
crossing between two ground states as a function of the 
exchange coupling J. The level crossing occurs in a sin- 
glet superconductor between states where the impurity 
spin is either unscreened or partially screened. We note 
here that the level crossing between two ground-state en- 
ergies does not occur at the point where the excitation 
level becomes a zero-energy excitation, since the free en- 
ergy of the superconductor is determined by taking into 
account continuum states as well as any localized intra- 
gap states in the self-consistent solution. 

We address the above problem at zero temperature 
by using the mean-field approximation both within the 
T matrix formulation and the self-consistent approach 
which takes into account a local gap-function relaxation. 
We also consider a local Coulomb interaction U which 
breaks particle-hole symmetry leading to an asymmetric 
spectral density for the impurity-induced quasiparticle 
states. Figure 1 illustrates the local effect of a magnetic 
moment on the low-energy spectral density in an s-wave 
superconductor. Since we limit our considerations to a 
classical spin, S ^> I, the impurity moment cannot be 
screened completely by the quasiparticles. Our main re- 
sults are as follows, (i) We show that the gross features of 
the impurity-induced quasiparticle states and their spec- 
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FIG. 1. The spin-unpolarized spectral density A(r = 0,lo) 
in an s-wave superconductor due to a magnetic moment lo- 
cated at r = as predicted by the T matrix approach. The 
continuum contribution is given in units of Nf and the height 
of the intragap bound-state peaks, located at fio = ±0.43Ao, 
denotes their integrated spectral weights in units of 7V_pAo. 
Here, Nf is the density of states at the Fermi energy in the 
normal state and Ao is superconducting energy gap. The 
magnitude of the magnetic moment is ttNfw = 0.7 and the 
local Coulomb interaction is ttNfU = 0.4. 

tral properties in s- and rf-wave superconductors can be 
qualitatively understood within the nonself-consistent T 
matrix formalism, (ii) For an s-wave superconductor, we 
find that, above the critical moment w c , the order param- 
eter changes sign at the impurity site and that it vanishes 
as w —> oo [13]. In the self-consistent mean-field approxi- 
mation, we find that order-parameter relaxation shifts w c 
downwards and that the energy of the impurity-induced 
bound state does not reach zero before a first-order tran- 
sition between the two ground states occurs, (Hi) In 
contrast, a rf-wave superconductor has no quantum tran- 
sition for any value of the magnetic moment when its 
quasiparticle spectrum in the normal state has particle- 
hole symmetry. The absence of the transition follows 
from the behavior of the impurity-induced quasiparticle 
states which are pinned at the chemical potential for an 
arbitrarily large magnetic moment. However, if particle- 
hole symmetry is broken or if the pairing state acquires a 
small s-wave component, the transition is again possible 
for a large enough moment. The rf-wave order parameter 
does not appear to change sign at the impurity site. The 
impurity moment induces two virtual-bound states which 
have four-fold symmetry and extend along the nodal di- 
rections of the energy gap. (iv) Finally, a mapping to an 
effective theory is derived which allows the formation of 
the impurity band at finite impurity concentrations to be 
explored. 

The Kondo effect in an s-wave superconductor inter- 
acting with a quantum impurity spin was considered in 
a recent work [15] which found a similar type of a quan- 
tum transition from the spin-doublet state (unscreened 
impurity spin) to the spin-singlet state (screened impu- 
rity spin). For the antiferromagnetic coupling (J > 0), 
the two states cross at T/r/Ao = 0.3, whereas, for the 



ferromagnetic coupling (J < 0), the coupling constant 
flows to the weak coupling limit and the level crossing 
is absent. The situation resembles the antiferromagnetic 
Kondo problem in that the moment is screened. How- 
ever, the crossover is replaced here by a sharp transition 
[15]- 

The physical picture of the quantum transition fol- 
lows from the behavior of the impurity-induced bound 
state. This transition is a consequence of the instability 
of the spin-unpolarized ground state, because, for a large 
enough w, the energy of the impurity-induced quasipar- 
ticle excitation would fall below the chemical potential. 
For a weak coupling w < w c , the ground state of the 
superconductor is a paired state of time-reversed single- 
particle states in the presence of the impurity scattering 
[16], |$o)»<» c ~ |$o), where 

|$o> = n( u .- + ^JT^-.-i)l°)' w 

and the bound state at the energy Qo is an intragap ex- 
cited state, \$-ii)w<w c ~ = tIijJ^o), where 

^-n} = ^UiU( u ' + v ^Ui)\°}- (2) 

i>l 

The quasiparticle operators [17] are 7^ = u^if — I'iV'lij 
and tIj'i = WiV'iij +Viipii', —i indicates the time-reversal 
conjugate of i and uf + vf = 1. As we find below, 
the time-reversal conjugate of this state, \<!>i-\} w < Wc ~ 
Tifl^o); does not appear in the superconducting energy 
gap. The processes of adding an electron with spin down 
to state i = — 1 and removing an electron with spin up 
from state i = 1 in the ground state \^o)w<w c lead 
to the same excited state: V'lijJ^o) = wil^-il) an d 
i/>i||$o) = Thus, in a tunneling experiment, 

the bound state is responsible through these processes for 
a finite differential conductance proportional to uf and vf 
at energies u> = ±Slo- Figure 2 elucidates schematically 
the mixing of particle and hole degrees of freedom in a 
quasiparticle state using a potential scatterer as an exam- 
ple. Because the particle and hole components are part 
of the same state, they appear symmetrically relative to 
the chemical potential. As w increases, Slo approaches 
the chemical potential and, at the critical coupling w c , it 
becomes a zero-energy state. At this point, the ground 
state of the superconductor becomes unstable. The new 
ground state of the superconductor becomes a paired con- 
densate for all but the impurity bound state with a fi- 
nite single-particle amplitude at the impurity -bound state. 

Thus, it has the form \^o)w>w c ~ | 14.) • This unpaired 

ground state is a result of the competition between the 
pairing-condensation energy and the magnetic interac- 
tion. For a strong magnetic interaction, the gain in the 
magnetic energy, due to the induced spin, dominates the 
condensation energy. The (variational) forms (1) and 
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(2) allow a varying number of doubly excited quasiparti- 
cles of the pure superconductor created by the operators 
t/|T1j'I from |$o)- These lead to a local suppression of 
the gap function A(r) and, in particular, for w > w c , to 
a negative value of A(r) at the impurity site (see below 
and Ref. [14]). The spin density is spread over a volume 
of linear dimension of the order of the coherence length. 

The plan of the paper is as follows. In Section II, the 
general formalism is introduced. Sections III and IV ex- 
plore the consequences of a magnetic moment in s- and 
d-wave superconductors. Both uniform and non-uniform 
order parameters are considered within the T matrix and 
self-consistent mean-field approximations. In Section V, 
the generalization to many impurities is outlined. Sec- 
tion VI contains a discussion of the number of electronic 
degrees of freedom associated with the impurity-induced 
quasiparticle states. In Appendix A, the interplay be- 
tween impurity scattering and pairing symmetry is fur- 
ther elaborated from the quasiparticle perspective. In 
Appendix B, the T matrix is generalized to finite impu- 
rity concentrations. 

II. FORMALISM 

Our starting point in describing a localized magnetic 
moment in a superconductor is the lattice formulation of 
electrons hopping between nearest-neighbor sites and in- 
teracting via an effective two-particle interaction. The 
pairing interaction is assumed to be weak or at most 
moderate so that the mean-field approximation gives a 
qualitatively reliable description of the superconducting 
ground state and the low-energy excitations. Specifically, 
consider a two-dimensional lattice model where a local- 
ized magnetic moment is created by a classical spin S at 
r = 0. The model is defined by the effective Hamiltonian 
H = Ho + -ffimp, where Ho describes a BCS superconduc- 
tor [18] and if; mp is the contribution due to the magnetic 
moment. In the mean-field approximation, 

(Rr) R 

-£A(R,r)tft(R + r)fi*(R), (3) 

Rr 

where (Rr) denotes nearest-neighbor sites separated by 
r, W is the half bandwidth (on a square lattice), \i is the 
chemical potential, and A(R, r) is the superconducting 

gap function. The operator \?(r) = [t/ , t( r ) V'{( r )] T i s a 
two-component Gor'kov-Nambu spinor, f a (a = 1,2,3) 
are the Pauli matrices for particle-hole degrees of free- 
dom, and fo is the unit matrix. Given that the pairing of 
electrons occurs in the spin-singlet state, the supercon- 
ducting order parameter (amplitude) can be expressed in 
the form 




FIG. 2. A schematic evolution of the spectral density in a 
superconductor as a function of a local Coulomb interaction 
U. The contribution due to the impurity-induced resonance 
state is depicted by the shadowed peaks. For \U\ ~ Ao, the 
resonance has nearly equal weights in the electron and the 
hole components, whereas, for \U\ 3> Ao, it is mostly either 
electron- or holelike. 

F(R, r) = \ J2(ih)*v(ipAR + r)^(R)). (4) 

av 

The relation between the order parameter and the gap 
function is given by the equation 

A(R, r) = —v(r)F(R,r). (5) 

Here, v(r) is the strength of the two-particle interaction, 
which is assumed to be instantaneous in time. Thus, the 
energy cutoff in the gap equation is set by the bandwidth. 
In a translationally invariant system, the BCS Hamilto- 
nian reduces to 

#o = ^'4(ekT3-A k f 1 )1' k . (6) 

k 

where A k = £ P A(R, r)e" ik r and * k = (^ kT ^l ki ) T . 
The fermion operators in real and momentum space 
are related by the unitary transformation, i/v( r ) = 
jy-1/2 t/> k(7 e 8k * r , where N is the number of sites in 
the system. For a square lattice with the nearest- 
neighbor hopping, the single-particle energy relative to 
the chemical potential in the normal state is e k = 
— -J W(cos k x a + cos k y a) — fi; a is the lattice spacing. 
At half filling, \i = 0, the effective Hamiltonian with 
nonzero on-site order-parameter is symmetric under the 
particle-hole transformation generated by the operator 
T\- *k —> * k = n*Q-k, where Q = (ir/a, ir/a). 



3 



The nature of the superconducting condensate can be 
established by studying the order parameter F which 
incorporates the pairing correlations and gives an es- 
timate for the size of the Cooper pairs [18]. (i) An 
on-site attraction favors s-wave pairing characterized by 
F(TL, r = 0) 7^ and A k = Ao. As an example, consider 
a uniform superconductor close to half filling (jjl ~ 0) so 
that the Fermi surface is square. Let n be the normal to 
the two sheets of the Fermi surface which are most nearly 
normal to r and rj_ = n • r > 0. Thus, rj_ > |r|||, where 
ry = (e3 x n) • r is defined. For rj_ a, 

F(R,r) = ^ Ko(r±/£, s ±)cos(k F _ L r_ L )6 ril o, (7) 

where A'o is the modified Bessel function of the zeroth 
order, £ s j_ = a(W/ Ao\/2) is the coherence length, and 
kp j_ ~ ir/a\/2. In the direction perpendicular to the 
Fermi surface, the extent of the pairing correlations is 
determined by the coherence length £ s j_, whereas, in 
the parallel direction, their range is vanishingly short, 
£ s || ~ 0. (n) For a strong repulsive on-site interac- 
tion and a nearest-neighbor attraction, pairing with d- 
wave symmetry may result [19]: F(TL, r = 0) = and 
At = 2Ad(cos k x a — cos k y a). Close to half filling and for 
rj_ a, one can again estimate [20] 



(ry/a) 



4tt \uJ [(m/«) 2 + K/^) 2 ] 3/2 



sin(fc FJ _rj_), 
(8) 



where £j_ = a(W/4:Ad) is the characteristic coherence 
length perpendicular to the Fermi surface, Ao ~ 4A<j 
is the maximum energy gap at the Fermi surface, and 
r) = sgn(r x r y ). This result is most useful when there is a 
clear separation between the length scales a and £j_ . For 
a given r j_ , its maximum value and the position scale as 
a ^J_/ r i an d r\\ = a(r j_ / '£ . The existence of nodes 
in the energy gap has resulted in the power-law decay of 
the pairing correlations in all directions at large distances 
instead of the exponential one as found in the s-wave 
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symmetry as expected. 



The interaction between the conduction electrons in 
the superconductor and the impurity spin is given by the 
Hamiltonian 



so that the z-axis points in the direction of w E ^^S. 
Obviously, the exchange interaction preserves particle- 
hole symmetry, whereas the Coulomb interaction breaks 
it. 

We will argue, after comparing the uniform order- 
parameter solution with the self-consistently determined 
one, that the difference between these two approaches 
is important only close to the impurity site. The as- 
sumption that the spatial variation of the superconduct- 
ing order parameter can be neglected leads to a con- 
siderable simplification of the problem, because the im- 
purity interaction (9) is limited to one site. Scattering 
of quasiparticles from the impurity moment is described 
by a T matrix, T(r,r';w) = T(u>)6 r o6 r io, which is lo- 
cal and whose Fourier transform therefore is independent 
of wave vectors. In the Gor'kov-Nambu presentation, 
the single-particle Green's function is G(r,r';t — t') = 
—i{T i S(r,t) i S^(r',t')). In the presence of an impurity, 

G(r, r'; u) = G (0) (r - r', u) + G ( °Hr, u)f(u)G^°\-T',u), 

(10) 

where both G^°\r,uj) = TV" 1 £ k G^(k, w)e ik r and 
T(u>) are matrices in two dimensional Gor'kov-Nambu 
space. The single-particle Green's function of a clean 
system is G^°\k,ui) = (uif - e k f 3 + A k fi + i0 + ) _1 . The 
T matrix can be obtained from the Lippmann-Schwinger 
equation in the form of 



f(w) = [T/- 1 -G(°)(r = 0,w)]- 



(11) 



These equations allow a complete solution of the problem 
as long as order-parameter relaxation can be neglected. 



III. MAGNETIC MOMENT IN AN S-WAVE 
SUPERCONDUCTOR 

Conventional pairing of s-wave symmetry is described 
by the gap function A(r) EE A(r, 0). It may be modeled 
by a system with an attractive on-site interaction v(r) = 
v(0)6 r0 , v(0) < 0. 



A. Uniform order parameter 



Hi 



imp 



JS • s(0) + Un(0), 



(9) 



N 



kk' 



where V = wfo + Ufs, J is the exchange coupling, 
U is the local Coulomb interaction, and the operators 

s ( r ) = \ Y,ov ^( r )^fi( r ) and n ( r ) = E CT V\Ur)Vv(r) 
are the conduction electron spin and number densities. 
Because of spin-rotational symmetry, the coordinate sys- 
tem for the spin degrees of freedom can always be chosen 



First, consider a superconductor at half filling with a 
spatially uniform s-wave gap function, A(r) = Ao. As 
shown by Yu [1], Shiba [2], and Rusinov [3], the magnetic 
impurity gives rise to a bound state inside the supercon- 
ducting energy gap. Their solution suggests that there 
exists a singular point, w = w c (U), that separates two 
distinct ground states of the system. The nature of the 
bound state can be found by computing the spectral den- 
sity per site 
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FIG. 3. (a) The bound-state energy fio, (b) the spectral 
weight Z^\o) (in units of AfAo), (c) the spin polariza- 
tion (s z (r = 0)} ((ra(r = 0)} = 1), and (d) the gap func- 
tion A(r = 0) at the impurity site r = as a function of w 
in the s-wave superconductor for U = 0. Lines denote the 
T matrix results for the uniform order parameter and sym- 
bols the self-consistent mean-field results on a square lattice 
at half filling. The quantities of the impurity-induced intra- 
gap quasiparticle state belonging to the branch w < w c are 
denoted by solid lines and solid symbols, whereas those ones 
belonging to the branch w > w c are marked by dashed lines 
and open symbols. Note that Z ( a \o) = Z ( ^j (0) because of 
particle-hole symmetry. 

4(r,u) = --ImG„(r,r;w). (12) 

7T 

For |w| < Ao, this can be expressed as A a (r,u>) = 
Z ( +\v)6{uj - fi ) + Z ( -\v)6{uj + fi ), where Z ( ±\v) 
describes the particle (hole) weight at positive (nega- 
tive) energy with spin orientation a (=T,I) and Slo is 
the energy of the bound state. Note that zi \r) = 
Z^J (r) whenever the normal-state energy spectrum ft 
has particle-hole symmetry and U = [21]. 

At half filling, /j, = 0, one can show that, for 
H < A , & °\r = 0,w) = -ttN(lo)(loto - 
A fi)/^A^ - ui 2 . For A < W, we obtain N(u>) = 
(2/tt 2 W) \og(4W/x/A 2 - lo 2 ). In the following, we will 
approximate N(u>) ~ N(0). The logarithmic dependence 
guarantees that this assumption is reasonable as long as 
|w| is not too close to Ao. Note that, for |w| <C W, 
N(0) equals to the normal-state density of states aver- 
aged over the energy interval ±Ao around /j, = 0. Thus, 
it is convenient to define the density of states per site 
in the normal state at the Fermi energy as Np = N(0). 
We can immediately solve the above equations for the 
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FIG. 4. (a) The bound-state energy fio and (b) the spectral 
weight ^^(O) (in units of 7V_pAo) of the impurity-induced 
intragap quasiparticle state as a function of U in the s-wave 
superconductor for ttNfw = 1.4. Lines denote the T ma- 
trix results for the uniform order parameter and symbols the 
self-consistent mean-field results on a square lattice at half 
filling. The quantities belonging to the branch w < w c (U) 
are denoted by the solid lines and solid symbols, whereas 
those ones belonging to the branch w > w c (U) are marked by 
dashed lines and open symbols. Note that Zi~\o) / Z ( ^J(0), 
because particle-hole symmetry is broken by the on-site 
Coulomb interaction U. We emphasize that while the spectral 
weights of the intragap peaks are asymmetric their energies 
relative to the chemical potential are not. 

bound-state energy and the spectral densities, (i) For 
w < w c (U) = ^(ttNf)- 2 + U 2 , the b ound-state energy 
Slo = ^< , where 

= r?A ^= C+C " ~ 1 (13) 
j{c\ + l)(c 2 _ + 1) 

c± = c w =p c u , and rj = sgn(c + c_). Here, the di- 
mensionless parameters Rrc Cyj — (w/ttN f )/(w 2 - U 2 ) 
and c u = (U/ttNf)/(w 2 — U 2 ). A straightforward cal- 
culation shows that z[ +) (0) = Z + SZ, Z ( f \o) = 
Z - 6Z, and ^ +) (0) = Z ( ~\{)) = 0, with Z = 
(A F Ao)2 7 r|cJ[(c 2 +c 2 ) + (c 2 -c 2 ) 2 ]/[(c?; + l)(c 2 _ + l)] 3 / 2 
and SZ = (N F A ) S gn(U)A7r\c u \cl /[(c\ + l)(c 2 _ + l)] 3 / 2 . 
Hence, the bound state is a local quasiparticle state with 
spin down [22]; i.e., the spin projection of the particle 
component is antiparallel and that of the hole compo- 
nent is parallel to the impurity moment w. Basic fea- 
tures of the low-energy spectral density are illustrated in 
Fig. 1. As m increases, the minimum energy to make 
a quasiparticle excitation from the condensate decreases, 
until, for w > w c (U), it would be negative. At this point, 
the superconducting condensate becomes thermodynam- 
ically unstable against spontaneous creation of a local 
quasiparticle excitation with spin down. As a result, the 
superconductor goes through a quantum transition from 
the spin-unpolarized to the spin-polarized state. Note 
that the average number of electrons remains the same. 
While, for U = 0, (n(r)} is a uniform function, (s z (r)} 
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becomes more peaked at r = as w increases. (For 
U ^ 0, also (n(r)) would be nonuniform.) However, the 
true nature of the transition is revealed by considering 
the total spin polarization, (s z ) = ^2 r (sz(^)} ■ At zero 
temperature, we obtain (s z ) = 0, for w < w c (U), and 
(s z ) = — 77, for w > w c (U). Because the ground state 
has qualitatively changed, so have the elementary exci- 
tations, and, in particular, the nature of the impurity- 
induced bound state, (ii) For w > w c (U), the energy of 
the bound state is given by Slo = ^>, with £1> = — £1<. 
Moreover, it is now associated with a quasiparticle with 
spin parallel to the impurity moment: ^"'"•'(0) = Z — SZ, 

Z ( ~\Q) = Z + SZ and z[ +) (Q) = Z ( ~\q) = 0. The 
bound state in both cases accounts for one electronic de- 
gree of freedom. 

Figure 3 summarizes the above results for U = 
and Figure 4 describes the effect of nonzero U . The 
Coulomb interaction effectively increases the critical cou- 
pling w c (U). Thus, for a large enough U, the doublet 
state is thermodynamically stable. It is also clear that 
the on-site Coulomb interaction breaks up the particle- 
hole symmetry. For example, given that the impurity mo- 
ment can be described classically and the envelop of the 
experimentally observed intragap spectral weight of the 
electronlike state Z^ + \y) is larger than that of the hole- 
like state Z(~\r), then either (a) U > and w < w c (U) 
or (b) U < and w > w c (U). 

Interestingly, one is tempted to interpret the recent ex- 
periment [11] on magnetic impurities as if the intragap 
spectral density decays on a length scale determined not 
by the coherence length but by a much shorter scale of the 
order of the Fermi length lp . This result appears para- 
doxical because a simple scaling estimate suggests that if 
the extent of the bound state were 0(£p), the energy of 
this state should be Tivp/ip ~ ep in a one-band picture 
and therefore could not be localized within the energy 
gap. One possible explanation of the experimental result 
is the power-law prefactor of the exponentially-decaying 
bound state. For example, consider a uniform gap func- 
tion, A(r) = Ao, and an isotropic energy band. In 
three dimensions, we find that Zf \v) oc {l F /r)' 2 e-' 2r l x °, 
where A = £o/\A - (fi /A ) 2 , £o = hv F /A , and 
£p = 2ir/kp. Thus, at the distance r = I0£p, Z^\r) 
has already decreased by a factor of ~ 100 from its max- 
imum value, which certainly should be enough to explain 
the experimental findings. Indeed, the spatial variation 
of the order parameter in the neighborhood of the impu- 
rity does not seem to play a major role in determining 
the behavior of Z^\r), even though the order parameter 
relaxes within a few Fermi lengths lp [24]. 



B. Non-uniform order parameter 

A qualitative understanding of the consequences of the 
order-parameter relaxation is obtained by considering a 
local suppression of the gap function at the impurity site, 
6 A = Ao — A(r = 0). Its effect can be accounted for by 
including the additional term 6A\?t(0)fi\Ij(0) in Hamil- 
tonian (9). Thus, V = wfo + Ufs + 6A?i . First, consider 
w = 0. The suppression of the order parameter provides 
a local attractive potential for quasiparticles which pro- 
duces two bound states in the superconducting energy 
gap at the energy 

fio = Ao^l-a 2 , (14) 

where a = C&//3. In addition, we have defined (3 2 = 
i(l - c 2 ) 2 + c 2 (J3 > 0), c 2 = c 2 - c 2 - 4, c w = 
(w/ttNf)/(w 2 - U 2 - 8A 2 ), c u = (U/ttNf)/(w 2 - U 2 - 
6A 2 ), and c A = (6A/TrN F )/(w 2 - U 2 - 6A 2 ). In con- 
trast to the bound state induced by the magnetic mo- 
ment, this state is two-fold degenerate by virtue of time- 
reversal symmetry. Their spin-quantum numbers are 
±7j. Note that the bound states become essentially in- 
distinguishable from the quasiparticle continuum (i.e., 
a 2 <C 1) in the limits of the weak order-parameter sup- 
pression, irNp 8A <C 1, and the strong potential scatterer, 
irNp \U\ ^> 1. Positive w breaks the degeneracy by low- 
ering the energy of the spin-down state and increasing 
the energy of the spin-up state: 

n /A = 1 -^fVT^± c -j^. (15) 

For w ^> 8A, the spectral density measured at the im- 
purity site r = is much larger for the former state at 
the energy u> < Ao than for the latter one whose spec- 
tral density is either transformed into the quasiparticle 
continuum or to the antibound state above the energy 
W. 

Next, compare the nonself-consistent T matrix predic- 
tions with the results obtained when the gap function 
is allowed to relax in the neighborhood of the impurity 
moment. Our numerical approach is based on the model 
defined by Eqs. (3) and (9) on a square lattice with the 
effective on-site electron-electron attraction of magnitude 
v(0) such that Ao/W = 0.05. The model, together with 
the mean-field equation, A(r) = — v(0)F(r, 0), is solved 
self-consistently (e.g., in the Bogoliubov - de Gennes 
scheme) [23]. As expected, the overall agreement between 
it and the T matrix calculation is good; see Figs. 3 and 
4. However, the additional variational degree of freedom 
shifts w c to a lower value, and Slo never becomes exactly 
zero. For w > w c , A(r = 0) is overscreened to a negative 
value and, ultimately, A(r = 0) —> 0, as w —> oo. The 
non-zero total spin polarization residing in the neighbor- 
hood of the impurity has lead to a nearly complete pair- 
breaking effect at this site. The charge density remains 
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FIG. 5. The spectral weight Z^ + \r) and Z^~\r) (in units 
of 7V_pAo) of the impurity-induced intragap quasiparticle state 
in the s-wave superconductor for (a) fi = 0, ttNfU = 0.4, 
and (b) fi = —W/2, U = 0. These results are computed 
self-consistently on a square lattice with the lattice spacing a, 
ttNfw = 0.7, and Ao/W = 0.05. The bound-state energies 
are (a) fi = 0.39A and (b) fi = 0.46A . 

uniform everywhere in the system for arbitrary w and 
U = 0. For w ~ w c (U), the results that assume a uni- 
form gap function everywhere deviate most clearly from 
those ones allowing the gap function relax locally, be- 
cause 8A/W and 8A/w have their largest values for the 
same values of w. 

At n = 0, the Fermi surface is a square producing a 
very anisotropic Fermi velocity. This anisotropy causes 
the impurity-induced bound state to have a shape which 
extends along the diagonal directions; see Fig. 5(a). A 
similar spatial form has also been observed in the context 
of relaxed quasiparticle excitations [25]. Notably, the de- 
cay of the spectral density due to the bound state is gov- 
erned by ^ n = hv a / Ao, where n is a normal to the Fermi 
surface at a given point and v a = 7} _1 (n • Vk£k)fc=fc F is 
the Fermi velocity at the same point. For a flat Fermi 
surface, zi^(r) oc 8 r ^e~ 2r - L l x± - , where ry and rj_ are 
the components of r parallel and perpendicular to the 
sections of the Fermi surface for which rj_ > |r|||. Here, 
Aj_ = <f s j_\/l — (Slo/Ao) 2 . This is in contrast to the 
isotropic r - 1 e ~ 2r / x o behavior found in two dimensions 
for a circular Fermi surface. The anisotropy in the Fermi 
velocity yields a similar structure in the gap function 
around the magnetic impurity [26]. The slower relaxation 
rate along the diagonal directions is clearly observable in 
Fig. 6(a). 

Moving away from half filling, the anisotropy in the 




FIG. 6. The self-consistent s-wave gap function A(r) 
around a magnetic moment for (a) fi = 0, ttNfU = 0.4, and 
(b) fi = —W/2, U = 0, obtained on a square lattice with the 
lattice spacing a. The other parameters are ttNfw = 0.7 and 
Ao/W = 0.05. 



Fermi velocity decreases, and the spatial form of the 
impurity-induced bound state and the gap function in 
the neighborhood of the impurity become increasingly 
isotropic, as shown by Figs. 5(b) and 6(b). Moreover, for 
H^zQ, the envelopes of the spectral weights ^"'"•'(r) and 

Z { ~\y) differ, b ecause the single-particle band structure 
of the normal state is no longer particle-hole symmet- 
ric. The sensitivity to the energy spectrum away from 
the vicinity of the chemical potential is a consequence 
of strong impurity scattering whose spatial extent is less 
than the lattice spacing and which efficiently mixes quasi- 
particle states with a wide wave- vector range around hp. 
A similar conclusion is reached in Ref. [14] where the 
effect of the band structure is discussed. 

The deformation of the order parameter in the vicin- 
ity of the magnetic moment creates an attractive poten- 
tial with a finite range. This feature may lead to addi- 
tional bound states which are no longer invariant under 
all symmetry operations of the square lattice (i.e., to 
higher "angular momentum" states) but whose binding 
energies are very small. Because these states have a node 
at the impurity site, they do not couple to the impurity 
moment and they are at least two-fold degenerate due 
to the spin degeneracy. In general, we find a two-fold 
degenerate state slightly below Ao that transforms like 
quadrupole (x 2 — y 2 ). However, away from half filling 
and for w < w c , there is some evidence for the lowest en- 
ergy, nonsymmetric quasiparticle state that is four-fold 
degenerate and transforms like dipole. 

IV. MAGNETIC MOMENT IN A D-WAVE 
SUPERCONDUCTOR 

Invariance of various perturbations under particle-hole 
transformation was significant in classifying the prop- 
erties of the impurity-induced quasiparticle states in s- 
wave superconductors. In a similar fashion, we may de- 
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scribe the same properties in rf-wave superconductors un- 
der particle-hole transformation. In rf-wave superconduc- 
tors, it is useful to introduce still another transformation, 
namely "charge conjugation": 'ft — ► = r^Q-k. 
It identifies the cases when the spin-unpolarized ground 
state of the superconductor interacting with a magnetic 
moment is unstable against spontaneous creation of a 
local spin polarization. As defined here, an s-wave su- 
perconductor is not invariant under charge conjugation. 



A. Invariance under charge conjugation 



£j( )(r = 0,w) at the chemical potential. This prevents 
S!± for crossing the chemical potential and ultimately 
pins it to zero energy as w —> oo with U fixed. 

In the neighborhood of the impurity, the virtual-bound 
state acquires the characteristic form of the rf-wave en- 
ergy gap with tails extending along the diagonal direc- 
tions. The spatial form of this state can be obtained 
by computing the spectral density, A a (r,u>) = Ao(u>) + 
8A a {v,u>), where Ao(u>) = — 7r _1 Im Go(lo) is the uniform 
spectral density of a clean superconductor and bA a (v,uj) 
is the impurity-induced contribution. For |w| <C Ao, 
Ao(u>) = Np |w|/Aq. The impurity-induced contribution 



First, we consider a quasi-two-dimensional system 
whose energy spectrum is particle-hole symmetric in the 
normal state and whose superconducting energy gap has 
d x 2_ y 2 symmetry at the Fermi energy: At = Aocos2<,s>. 
We also assume that the gap function is uniform every- 
where and that the Fermi surface in the normal state is 
circular. Since the system possesses charge-conjugation 
symmetry, T2&°' > (r = 0,w)f2 = &°\r = 0,w), it is im- 
mediately clear that &°\r = 0,w) = Go(w)to. At the 
impurity site, the spectral density has a particularly sim- 
ple form: 



6A a (v,w) = --Im[G ( °\T,u)T!f\u)G ( +\-T,u) 



A a (r = 0,u) 



-— Im — 

7T 1 



av a G (oj) ' 



(16) 



where v± = w±U,a = + for spin up, and a = 
— for spin down. The poles of A a (r = 0,w) give 
information about the quasiparticle states. At low 
energies |w| <C Ao, one obtains the formula [5,6] 
Go(io) = -# F (2 W /A )[log(4A /|w|) + isgn( W )7r/2]. For 
l c ±l *C 1, where c± = (ttNfv±)~ 1 , the magnetic 
moment induces two quasiparticle states at energies 
S!± = -|c±Ao/log(8/7r|c±|) with inverse lifetimes T± = 
■||fi±|/log(8/7r|c±|). For U = 0, the two states are de- 
generate: = Their spin quantum numbers are 
s z = — i. With increasing U, the degeneracy is lifted as 
decreases and Sl_ increases. For U w, they appear 
as resonances in the spectral densities A^{v = 0,w) and 
j4|(r = 0,w) at u> = — £1+ (holelike) and u> = Sl_ (elec- 
tronlike), respectively. For U ~ w, there is no solution for 
Sl_ and no resonance structure develops in A^(r = 0,w). 
For U w, Sl_ is negative and approaches zero from 
below. Thus, for each spin orientation, the spectral den- 
sity A a (r = 0,w) has a holelike resonance at the energy 
ijj = — Sl<j. However, the quasiparticle energies £1± never 
cross the chemical potential for any value of w and U . 
This ensures that the superconducting ground state re- 
mains stable against spontaneous creation of quasipar- 
ticle excitations. Thus, in contrast to the s-wave case, 
a superconductor with rf-wave pairing does not exhibit a 
transition where the spin-quantum number of the ground 
state changes from zero to The difference follows from 
charge-conjugation symmetry which leads to vanishing 



+ G { °\y,w)T^\w)G^\-y,w)\, (17) 

where T ( ± \lo) = -l/[G±(w) - vvgl]. We use the 
notation in which G^\tc,^) = ^Tr t a &°\r, ui) and 
G a (u>) = G^a\? = 0,w) (a = 0,...,3). In addition, 
we have defined G^ 0) (r,w) = G^\r,u) ± G^\r,u) and 

G±(ui) = G±\r = 0,w). Note that particle-hole sym- 
metry implies the relation G±(u>) = Go(w). For U = 0, 
T+\ui) = T^\u>), and the spectral density is particle- 
hole symmetric: 6A^(r, —u>) = 6A^(r,u>). While the 
overall spatial dependence of the virtual-bound states 
generated by purely magnetic (U = 0) and nonmagnetic 
(w = 0) scattering is the same, the resonance factors 
lead to distinct states. The specific form of the 
spectral density is determined by the Green's function, 
G (0 \r,ui). Eq. (17) prod uces a qualitatively correct de- 
scription of the quasiparticle states as long as the energy 
gap is not comparable to the bandwidth, because then 
the relaxation of the gap function in the vicinity of the 
defect can be neglected as a small correction. Nonethe- 
less, the local suppression of the gap function A(r = 0) 
could easily be included. For u> = I f£l±, the spectral 
density is most strongly enhanced in the horizontal and 
vertical directions (the extrema directions of the rf-wave 
energy gap) at distances r <C £o (= hvp/Ao) from the 
impurity. Further away from the impurity (r ~ £o); the 
situation is reversed: the spectral density is enhanced in 
the neighborhood of the diagonal directions at the reso- 
nances. 

That the Green's function &°\r = 0,w) is propor- 
tional to fo has curious consequences which set rf-wave 
pairing apart from s-wave one. In the s-wave super- 
conductor, low-energy quasiparticles have both parti- 
cle and hole character due to the non-zero pairing field 
F(TL, r = 0). In contrast, rf-wave symmetry forces the 
pairing field to vanish locally, F(H, r = 0) = 0: the par- 
ticle and hole degrees of freedom are effectively decoupled 
at the impurity site. Strong magnetic scattering leads to 
two virtual bound states whose components at the impu- 
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FIG. 7. The self-consistent d-wave gap function [28] A(r) 
around a magnetic moment (ttNfw = 10 and {7 = 0) for 
(a) /t = and (b) /t = — IT/2 obtained on a square lattice 
with the lattice spacing a and the coherence length fj_ = 10a 
(4Ad/W r = 0.1). The minimum at r = is not shown. Ad is 
the spatially uniform d-wave gap function in the clean system. 
Note that the maximum energy gap at the Fermi surface is 
(a) A = 4A d and (b) A = 2A d . 

rity site are either electronlike with spin down or holelike 
with spin up; the other components vanish because their 
orbital character is the same as that of the order param- 
eter (i.e., they have a node at the impurity site). A non- 
zero U lifts the degeneracy without mixing the states. 
Note that, for an impurity whose non-magnetic scatter- 
ing strength (U > 0) is larger than the magnetic one, 
two virtual-bound states of hole character with quantum 
numbers s z = ±77 are obtained at the impurity site [6]. 

In conclusion, a rf-wave superconductor appears locally 
as a simple metal with a pseudogap: at short distances, 
the particle and hole excitations are decoupled and the 
energies of the virtual-bound states do not cross the 
chemical potential due to the level repulsion. Another 
treatment in terms of elementary quasiparticle excita- 
tions is outlined in Appendix A. 

B. Broken charge-conjugation symmetry 

First, assume a mixed s- and rf-wave superconductor so 
that the pairing state is described by a gap function At = 
Ao cos 2ip + A s , where a small s-wave component has 
been added to the rf-wave state, A s <C Ao. This type of 
admixture of two angular-momentum states is relevant in 
a d-wave superconductor where the tetragonal symmetry 
is broken by an orthorhombic distortion. It also breaks 
charge-conjugation symmetry. Under these conditions, 
G(°)(r = 0,w) = GoHfo + GiHfi; for M < A , G (u) 
is given above and G\{u>) ~ irNp A s /Ao. For purely 
magnetic scattering (U = 0), the energies of the virtual- 
bound states, S!± = -|2±Ao/log(8/7r|2±|), where c± = 
c w ± A s /Ao, are split by the s-wave component. Clearly, 
the state at u> = Sl_ crosses the chemical potential for 
c w < c* = A s /Aq. The critical coupling c* signifies a 




FIG. 8. The spectral density <5A CT (r, ±Qo) (in units of Nf) 
for (a) fi = and (b) fi = —W/2 in a two-dimensional d-wave 
superconductor as a function of position r = (x,y) around 
a classical magnetic moment (ttNfw = 10 and {7 = 0) lo- 
cated at r = 0; a is the lattice spacing. These results 
are computed self-consistently with fj_ = 10a. At half 
filling, the spectral density obeys particle-hole symmetry: 
SAf(r,—Qo) = <*>Aj(r, Q ). The energies of the shown vir- 
tual-bound states are (a) fio — 0.05Ao and (b) fio — 0.5Ao. 

quantum transition from the spin-unpolarized to a spin- 
polarized state of the superconductor. We remark that, 
for an s-wave superconductor, c* = 1. As A s — ► 0, the 
critical coupling c* — ► 0, and the quantum transition is 
shifted to higher values of the magnetic moment until, 
for A s = 0, the transition is completely suppressed. 

Second, if the particle-hole symmetry is broken in the 
normal state, the Green's function acquires a non-zero 
T3 component. For a large enough magnetic moment, it 
will again induce a transition where the spin-quantum 
number is changed from zero to -|. 

In general, both of the above features may be present 
simultaneously. In such a case, one must replace c± by 
c w ± \/ q\ + (cu + 93) 2 , where g a = (7rA F )" 1 G„(fi) (a = 
1,3) and w > \U\. Similarly, c* is replaced by c w for 
which c_ =0. As a consequence, the critical coupling 
becomes a function of the chemical potential and the s- 
wave component: c* = c H ,( / u, A s ). 

C. Non-uniform order parameter 

The Fermi- velocity distribution affects the precise form 
of the d-wave gap function [27] in the neighborhood of the 
magnetic moment; see Fig. 7. As in the s-wave super- 
conductor, the sensitivity to the Fermi-surface geometry 
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originates from the mode composition of the gap func- 
tion: the electronic degrees of freedom within the energy 
proportional to Ao about the chemical potential have the 
largest weight in the gap function. Close to half filling, 
the gap function has four troughs extending along the 
diagonal directions, whereas, away from half filling, the 
four-fold symmetry of the energy gap is more crucial in 
determining the spatial relaxation of the gap function. 
That the quasiparticle states are concentrated along the 
horizontal and vertical directions at energies for which 
the density of states is the highest (w ~ Ao) explains the 
latter finding. In the strong-coupling limit, some modifi- 
cations are expected [29]. 

A strong magnetic moment creates virtual-bound 
states whose spatial dependence has distinctive features 
due to d-wave pairing symmetry. For example, in the 
neighborhood of the defect, both the particle- and hole 
components of the spectral density remain non-zero along 
the diagonals, although the maximum value at a given 
distance is located off-diagonally; see Fig. 8. Close to 
half filling, the Fermi-surface geometry has a strong fo- 
cusing effect on the spatial variation of the virtual-bound 
states, and their lobes along the vertical directions spread 
only slowly. At half filling (jjl = 0), the lattice Green's 
functions G^\v,lj) and G^\r,ui) vanish identically at 
the sites for which (x + y)/a is an even integer. Because, 
for |w| <C Ao, G^(r,u>) is very small, the spectral den- 
sity is negligible along these lines. Note that, in general, 
Gi*\r,u>) does not vanish at the same lattice sites for 
s- and d-wave pairing, explaining the difference between 
Figs. 5(a) and 8(a). The spatial variation of the virtual- 
bound states is accurately given by Eq. (17) even though 
it assumes a uniform gap function as one can verify in 
this case. For rj_ a and U = 0, one can estimate 

8A a (r,u ~ 0) oc sin 2 (fc F _ L r_ L )/[(r||/a) 2 + (r ± /Uf], 

(18) 

where rj_ > ry (by definition). At half filling, kp±r± = 
ir(x + y)/2a. Away from half filling (ji = — W/2), the 
spectral density tends to vary radially demonstrating 
more clearly the effect of the k-space structure of the 
superconducting energy gap. The example shown in 
Fig. 8(b) is computed in the spin-polarized ground state 
which is obtained when the magnetic moment is large 
enough and the superconductor is not half filled. Neither 
does the spectral density obey particle-hole symmetry in 
this case. In general, quasiparticle excitations are quali- 
tatively described by Eq. (17) which provides a straight- 
forward method for determining the spectral density and 
the resonance energies. For example, the resonance en- 
ergies are obtained from the equation, G±(S!) = l/v±. 
The two-fold degeneracy of the virtual bound state is 
lifted both by the local Coulomb interaction and by the 
band structure which is asymmetric relative to the chem- 
ical potential. Because the lattice Green's functions are 



0.06- 6.0- 




(a) (b) 



FIG. 9. The electron spin density (s z (r)} in a 
two-dimensional d-wave superconductor as a function of 
position r = (x,y) around a classical magnetic moment 
(ttNfw = 10 and {7 = 0) located at r = 0; a is the lat- 
tice spacing. The chemical potential is (a) fi = and (b) 
/i = —W/2. The peak (s z (r = 0)} ~ — i in both cases is 
cut off in order to illustrate the finer details. These results 
are computed self-consistently with fj_ = 10a. Note that (a) 
(s z ) = and (b) (s z ) = -§. 

fast to compute numerically, we have been able study in 
detail the evolution of the spectral density as a function 
of the chemical potential and the impurity potential (not 
shown). 

In spite of the possibility of the quantum transition, 
A(r = 0) does not change sign for any value of mag- 
netic moment w. While, for c* < c w <C 1, the electron 
spin density is nearly completely polarized at the impu- 
rity site, {s z (r = 0)) ~ — -|, there is a compensating 
spin-density cloud of opposite sign in the neighborhood 
of the magnetic moment that leads to the vanishing spin 
polarization. Figure 9 shows the spatial variation of the 
electron spin density with a given magnetic moment both 
at and away from half filling, illustrating the two cases 
c w > c*(n) (ji = 0) and c w < c H ,( / u) (ji ^ 0). The mag- 
netic moment perturbs electronic degrees of freedom on 
the energy scale of w which must be large compared to 
the superconducting energy gap Ao to have pronounced 
virtual-bound states. At half filling, the behavior of the 
spin density is dominated by the anisotropy in the Fermi 
velocity, disguising the effect of pairing symmetry. Away 
from half filling, the virtual-bound states contribute at 
distances r > £o; while at shorter distances the non-zero 
spin density is mostly coming from quasiparticle states 
residing outside the energy-gap region. 

These considerations offer intriguing prospects regard- 
ing high-temperature superconductors. In addition to 
possible d-wave features, their low carrier densities ex- 
acerbate the role of the highly anisotropic Fermi-surface 
geometry that is inevitably mixed in the spatial depen- 
dence of the gap function and the quasiparticle states in 
the vicinity of impurities. These properties should be ac- 
cessible to many experimental probes, such as scanning- 
tunneling microscopy, etc. In particular, it might now 
be possible to determine experimentally the coherence 
length in short-coherence-length superconductors by di- 
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rectly measuring the spectral density in the vicinity of 
defects. In low-temperature superconductors, this is not 
possible because the power-law prefactors kill the signal 
before the exponential cutoff set by the coherence length 
can be observed. Even though the coherence length in 
these materials is typically short (£o ~ 20 A), it still is 
long enough that our mean-field results should be rele- 
vant in describing their properties qualitatively [30]. 

V. MANY MAGNETIC IMPURITIES 

As an example, consider strongly-scattering magnetic 
impurities in an s-wave superconductor. In the dilute 
limit of impurities, only nearest-neighbor impurity states 
interact due to their exponentially decreasing overlaps. 
Appendix B describes a mapping of a system of dilute 
impurities onto an effective model. Using known re- 
sults for such a model, we conclude that all quasipar- 
ticle states in the impurity band in two dimensions are 
localized. In three dimensions, there is a critical density 
of impurities, below which all states are localized and 
above which a mobility edge appears. A similar prob- 
lem of the impurity-induced quasiparticle states states in 
a two-dimensional rf-wave superconductor has also been 
considered [31]. 

Although the assumption that the impurities do not 
interact leads to an equal density of states at every impu- 
rity site, the fact that the impurity-induced quasiparticle 
states may have large overlaps modifies this conclusion. 
Disorder with strongly interacting impurities yields large 
impurity- and energy-dependent variations in the quasi- 
particle wave functions. These variations may, in the 
worst case, obscure the detection of bound states — for 
example, by scanning tunneling microscopy — or facili- 
tate the detection if the probe is close to impurity sites 
that have very large relative weights in the wave function. 

VI. FINAL REMARKS 

There is a qualitative difference between s- and d- 
wave superconductors regarding the quasiparticle states 
induced by magnetic moments. This difference arises 
from the local properties of the superconducting order 
parameter F(H, r) which is non-zero for s-wave pairing 
but vanishes for rf-wave pairing at r = 0. Thus, in the 
s-wave superconductor, any local probe sees particle and 
hole degrees of freedom that are coupled by the pairing 
field. This feature suggests that the two intragap peaks 
in the spectral density should be interpreted as reflections 
of the same quasiparticle state. In contrast, because the 
pairing field in the rf-wave superconductor vanishes at 
short distances, a local perturbation leads to two sepa- 
rate quasiparticle states which are either purely particle- 



or purely holelike at the impurity site. One can distin- 
guish the two cases by their response to a local Coulomb 
interaction U n(0). The properties of the quasiparticle ex- 
citations are affected by U because their charges are not 
the same. Using the ground-state charge as a reference 
point, define the charge operator Q as 

Q = $>(r)-($o|n(r)|$o>]. (19) 

r 

Clearly, ($_ijJQ|$_i^) = u\ — v\. If the particle and hole 
features belong to two distinct states, the Coulomb inter- 
action U will shift their energies differently because their 
charges differ. In contrast, if they are manifestations of 
the same state, the Coulomb interaction could not affect 
them independently, and the spectral density will have 
two intragap peaks that are symmetrically located be- 
low and above the chemical potential. Indeed, we always 
find that, in an s-wave superconductor, the spectral den- 
sity contains two peaks at energies ±Slo regardless of the 
local Coulomb interaction whereas, in a rf-wave supercon- 
ductor, the Coulomb interaction shifts the quasiparticle 
peaks asymmetrically relative to the chemical potential. 

We conclude by noting that the special nature of the 
induced quasiparticle states has interesting implications 
for optical absorption in an s-wave superconductor con- 
taining a local moment. Because the bound state inside 
the energy gap exists only for a given spin direction (its 
time-reversed conjugate overlaps with the quasiparticle 
continuum), no sharp absorption feature is found at 2£lo 
but the optical absorption begins at the energy Ao + Slo- 

ACKNOWLEDGMENTS 

We would like to thank J. Byers and M. Flatte 
for interesting discussions and, in particular, for help- 
ful comments on the Fridel oscillations in the normal 
state, A. Yazdani for discussions on experimental issues, 
C. Kallin for introducing M.I.S. to the Bogoliubov - de 
Gennes technique, and A.J. Berlinsky for useful com- 
ments on the manuscript. This work was supported 
in part by Natural Sciences and Engineering Research 
Council of Canada, the Ontario Center for Materials 
Research (M.I.S.), by the U.S. Department of Energy 
(A.V.B.), by the NSF Grant No. DMR-9629987 (J.R.S.), 
and by the Many Body Theory Program at Los Alamos. 

APPENDIX A 

The reasons for different behavior of various types of 
impurities in superconductors, as well as the effect of 
pairing symmetry, can be illuminated by rewriting the 
problem in terms of the elementary quasiparticle exci- 
tations of the clean superconductor, described by the 
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Hamiltonian, Eq. (6). This Hamiltonian is diagonalized 
by the Bogoliubov-Valatin transformation, \? k — ► r k = 
5k*k, where T k = (7 kT 7l k |) T and 

5 k =f C ° S 2 k - Sm M. (Al) 
y sin t> k cos t> k J ' 

In this parametrization, i/ k = cos # k and i> k = sin^. 
Choosing sin26' k = A^/E^ gives 

ff = X)^ r k^r k , (A2) 

k 

where E^ = ^/e k + A k . In terms of these excitations, 
the Hamiltonian for the (^-function impurity is if; mp = 

J2a H(a) > with 

H(a) = ^J2 T i v ^-- ( A3 ) 

kk' 

The three types of scattering interactions (a = w,u,A) 
are: (i) magnetic, 

V&) = w[cos(9 k - k ,)f o - sin(# k - k ,)if 2 ]; (A4) 
(n) potential, 

^kk< = U[cos(9 k + k ,)f 3 + sin(0 k + k ,)fi]; (A5) 
and (Hi) local order-parameter suppression, 

^ = -<5A[sin(# k + k ,)f 3 - cos(0 k + k ,)fi], (A6) 

where 6A = Ao — A(r = 0) > 0. These expressions form 
the basis for the following discussion of when to expect 
a bound state in the superconducting energy gap. 

1. s-wave superconductor 

First note that, for k ~ hp, k ~ ^-sgn(A k ). For an 
s-wave superconductor at energies close to the chemi- 
cal potential, the magnetic interaction term is diagonal: 
V^} ~ wfo- Thus, the excitation energy is locally de- 
creased for spin-down excitations and increased for spin- 
up excitations (recall that the effective w changes sign 
for opposite spins of quasiparticles). As a consequence, 
a bound state appears only for the spin-down quantum 
number [32]. For large enough w, the bound-state energy 
crosses the chemical potential. 

The potential scattering term is purely off-diagonal: 
Vtf$ ~ Uf\ ; it causes level repulsion regardless of the sign 
of U . Therefore, potential scattering does not introduce 
intragap states. 

Local suppression of the order parameter leads to an 
attractive diagonal potential for both type of quasipar- 
ticles: V^ kk V ~ — 6Af 3 . This produces two degenerate 
bound states in the energy gap, as expected from time- 
reversal symmetry. 



2. d-wave superconductor 

In contrast to the s-wave case, the alternating sign of 
the d-wave order parameter yields quasiparticle scatter- 
ing terms that are predominately neither diagonal nor 
off-diagonal. While a large enough potential will lead 
to pronounced resonance states, they do not cross the 
chemical potential when the normal state quasiparticle 
spectrum has particle-hole symmetry and the gap func- 
tion average over the Fermi surface vanishes. 

Finally, the appearance of degeneracies in the spec- 
trum of virtual bound states can be understood in terms 
of the symmetry properties of the effective Hamilto- 
nian, H = Ho + -ffimp- In general, a rf-wave supercon- 
ductor is invariant under the transformation in which 
^ic — ► \? k = f 3 \?i{ k , where R denotes a rotation through 
an angle -| about the impurity site. When the system is 
half filled and U = 0, the symmetry group of H includes 
the charge-conjugation transformation. Because such a 
symmetry group is non-Abelian, it also has higher than 
one-dimensional irreducible representations. In particu- 
lar, the impurity-induced virtual bound states belong to 
a two-dimensional irreducible representation; i.e., they 
are two-fold degenerate. The degeneracy is lifted by a 
nonzero U, by an asymmetric band structure, or by a 
gap function which has a nonzero s-wave component. For 
w = 0, the quasiparticle states are at least two-fold de- 
generate because of time-reversal symmetry. 

APPENDIX B 

In this Appendix, a mapping to an effective theory 
is presented, from which the consequences of a finite 
density of impurities can be explored. To formulate 
the problem of many magnetic impurities, we utilize 
the four-dimensional Gor'kov-Nambu representation, be- 
cause the magnetic moments are allowed to be randomly 
oriented. In this representation, the effective Hamilto- 
nian is H = Ho + -ffimp , where 

#0 = ^*k( e kT3 + A k f 2 ^ 2 )* k , (Bl) 

k 

and \? k = (t/'kt V'-kj. V'if i J ^-ki) T is the Gor'kov-Nambu 
spinor, f a and a a (a = 1, 2, 3) are the Pauli matrices for 
particle-hole and spin degrees of freedom, and fo and &o 
are the unit matrices. We further assume that randomly 
distributed impurities scatter by (^-function interactions, 

H lmp = J2^(r)V(r)^(r) (B2) 

r 

where 

V(r) = J2v n 6rr n - (B3) 

n 
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Here, v n = w n -e+U n fs is the interaction strength of the 
nth impurity moment located at r n . Following the usual 
convention, we have defined e = (72^2 + T'3(o'iei + 17363). 
In the absence of the impurities, the single-particle 
Green's function is 



zk r 



TV to - e k f 3 - A k f 2 ^ 2 + i0 + ' 

k 



(B4) 



Generally, the total Green's function in the presence 
of impurities is given by 

G(x,x') = G (0) (x-x 1 ) + (B5) 
fdydy'G (0 Hx-y)T(y,y')G (0 Hy' -x'), 



superconductor. For example, for an s-wave supercon- 
ductor, Ao(u>) = N f\uj\ / \/ uj' 2 — Ag, when |w| > Ao, and 
zero otherwise. The second term, 



Im y {a 



|G(°)(i 



(Bll) 

uj)f mn {uj)G (0 \Y n -Y,uj)\a), 



incorporates the effect of impurities. We have adopted 
the notation, according to which (<t|G|<7} = G aa , etc. 
(Note that G is the particle component of G: G = Gn.) 
Because the impurities modeled by Eq. (B3) do not 
change the total number of states, Jdu> 8A a (r, u>) = 0. 
Eq. (Bll) can be written in a compact form, if we con- 
sider the spectral density only at a given impurity site 



where T is a T matrix. We use the notation in which 
all the matrices are given in four-dimensional Gor'kov- 
Nambu space, x = (r,t), and Jdx = ^2 r fdt. The T 
matrix is a solution of the Lippmann-Schwinger equation, 

f(y, y') = V(y)6(y - y') + ( dz V(y)G (0 \y - z)f(z, y'). 



(B6) 

For impurities with (^-function interactions, it is straight- 
forward to see that a T matrix of the form 



T(r,r';w) = ^ f mn (u>)6 rrm 6 r - rn , 



(B7) 



solves the Lippmann-Schwinger equation. Indeed, defin- 
ing 



(m|T(w)|n) = T m „(w), 
{m\G (0 \u 1 )\n) = G (0 \r m -r n ,u 1 ), 
(m|v _1 |n) = v'* ra „; 

the solution is 

T(w) = [V 1 - G(°)(w)]- 1 . 



(B8) 



(B9) 



This formula allows us to compute various physical prop- 
erties of the system in the case of a finite number of ran- 
dom impurities for s- or rf-wave superconductors. One 
such quantity of interest is the local spectral density. 
Eq. (B9) also implies a useful correspondence between 
the current problem and that of noninteracting quasipar- 
ticles moving on an infinite lattice with random on-site 
energies and hopping amplitudes (see below). 
Consider the local spectral density, 



1 



A a {v,uj) = Im(? CTCT (r,r;w), 



(BIO) 



and define A a (r,u>) = Ao(u>) + 8A a (r, u>), where Ao(u>) = 
— 7r _1 Im Gct°ct (r = 0,w) is the spectral density in a clean 



6A a (r n ,u) = --lm{na\G^\ui)T(ui)G^\ui)\na). 

IT 

(B12) 

The importance of Eq. (B12) becomes clear by a fol- 
lowing example. Assume that for a given u> = u„ there 
exists a state \<p a ) such that 



G(°)(o 



Wo 



0. 



(B13) 



Then, averaging over the random orientation of magnetic 
moments, the spin-unpolarized spectral density becomes 

8A{v n ,uj) = C n 'Y^ l {ncr\ip a ){ip a \n(T)C a {Lo), (B14) 



where G n = (w 2 + U%)/(w 2 — U 2 ) 2 and the Lorentzian 
function C a (u>) is given by C a (u>) = (27r) _1 / a w"/[(w — 
oj' u ) 2 + oj" ], u) a = u)' a + iu>", and /„ is the residue of 
the ath pole. This form is useful when /„ depends only 
weakly on u>. The state \<p a ) specifies the properties of 
quasiparticles at a given energy u> a . In particular, the 
matrix element |(n<7|<,s> a }| 2 incorporates the spatial in- 
formation of the impurity-induced states; for example, 
whether the quasiparticle state at the energy u> = u„ is 
localized or extended. 

The above result suggests a fruitful formulation for ex- 
ploring impurity-induced quasiparticle states in the su- 
perconducting energy gap. The central idea is to define 
the Hamiltonian, 



H 



7 y c m v m c m 



(B15) 



where c m is a four-component spinor and i mn = 
&°\r m — r„;w). Noting the similarity with Eq. (B13), 
we can immediately conclude that the zero-energy eigen- 
state of % also determines the properties of a quasipar- 
ticle state in the impurity band at a given energy u>. 
For example, in a three-dimensional s-wave supercon- 
ductor where magnetic impurities yield localized states 
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in the energy gap, the formation of the impurity band 
is mapped to the Hamiltonian % with short-range hop- 
ping amplitudes. Indeed, t mn oc j» _1 e _r ""»/ A , where 
r mn = |r m — r n |. In rf-wave superconductors, Eq. (B15) 
has been applied to study localization properties of quasi- 
particle states that are induced by unitary scatterers [31]. 



[1] L.Yu, Acta Physica Sinica 21, 75 (1965). 

[2] H.Shiba, Prog. Theor. Phys. 40, 435 (1968). 

[3] A.I. Rusinov, Sov. Phys. JETP Lett. 9, 85 (1969). 

[4] P.C.E. Stamp, Journ. of Magnetism and Magn. Mat. 63- 
64, 429 (1987). 

[5] A.V.Balatsky, M.I.Salkola, and A.Rosengren, Phys. Rev. 
B51, 15547 (1995). 

[6] M.I.Salkola, A.V.Balatsky, and D. J.Scalapino, Phys. 
Rev. Lett. 77, 1841 (1996). 

[7] Y.Onishi et al., Journ. Phys. Soc. Japan 65, 675 (1996). 

[8] L.Dumoulin, E.Guyon, and P.Nedellec, Phys. Rev. B16, 
1086 (1977); and references therein. 

[9] H.F.Hess, R.B.Robinson, R.C.Pynes, J. M. Varies, Jr., and 
J.V.Waszczak, Phys. Rev. Lett. 62, 214 (1989). 
[10] I.Maggio-Aprile, Ch.Renner, A.Erb, E.Walker, and 

0. Fischer, Phys. Rev. Lett. 75, 2754 (1995). 
[11] A. Yazdani, B.A.Jones, C.P.Lutz, M.F.Crommie, and 

D.M.Eigler (unpublished). 
[12] A.Sakurai, Progr. Theor. Phys. 44, 1472 (1970). 
[13] During the course of preparing this manuscript, we be- 
came aware of work by Flatte and Byers [14] also con- 
sidering magnetic impurities in s-wave superconductors. 
Our results are partially in accordance with theirs in 
that the structure of the s-wave gap function is distorted 
and it changes sign at the impurity site at the transition 
point. 

[14] M.E. Flatte and J. M. Byers (unpublished). 

[15] O. Sakai et al., Journ. Phys. Soc. Japan 62, 3181 (1993); 

and references therein. See also J.Zittarts and E.Miiller- 

Hartmann, Z. Phys. 232, 11 (1970). 
[16] P.W.Anderson, J. Phys. Chem. Solids 11, 26 (1959). 
[17] J.R.Schrieffer, Theory of Superconductivity (Benjamin- 

Cummings, Reading, 1983). 
[18] J.Bardeen, L.N.Cooper, and J.R.Schrieffer, Phys. Rev. 

108, 1175 (1957). 
[19] V.J.Emery, Phys. Rev. Lett. 58, 2794 (1987). 
[20] In this expression, a term of order C(a/{j_) 2 relative to 

the maximum value of F is omitted as a small correction; 

it is obtained by interchanging ry and rj_. 
[21] This relation follows from particle-hole symmetry: 

hG^Hr = 0,Lo)h = &°\r = 0,lo). 
[22] Because the spin projection along the impurity moment 

jdvs z (v) describes a conserved quantity, the quasiparti- 

cle states can be labeled according to the spin quantum 

number s z . 

[23] See, for example, P.G. de Gennes, Superconductivity of 
Metals and Alloys, (Addison- Wesley, New York, 1989). 

[24] R.Kummel, Phys. Rev. B6, 2617 (1972). 

[25] A.R.Bishop, P.S.Lomdahl, J.R.Schrieffer, and 
S.A.Trugman, Phys. Rev. Lett. 61, 2709 (1988). 

[26] A non-magnetic impurity has an analogous effect on the 



s-wave gap function; see T.Xiang and J.M.Wheatley, 
Phys. Rev. B51, 11721 (1995). 

[27] In our numerical mean-field calculations, a superconduct- 
ing ground state of d-wave pairing symmetry is realized 
for an effective two-particle interaction which is repulsive 
for electrons at the same site and attractive between the 
nearest-neighbor sites: v(r) = voS r o + v\ i r r„, where 
vo = — vi > and v labels the four nearest-neighbor 
lattice vectors on a square lattice. 

[28] The d-wave gap function A(r) is defined as the ampli- 
tude of the partial-wave component that transforms as 
d x 2_ y 2 under the symmetry operations of the square lat- 
tice: A(r) = j^2 v ( — l) v /\(v,v v ), where v sums over the 
four nearest-neighbor lattice vectors and r_„ = — r„. Be- 
cause the magnetic moment breaks up locally the crystal 
symmetry, also other partial waves acquire non-zero am- 
plitudes in the vicinity of the moment. 

[29] M.Franz, C.Kallin, and A. J.Berlinsky, Phys. Rev. B54, 
R6897 (1996). 

[30] In contrast to Ref. [29], which focuses on non-magnetic 
defects in strong-coupling d-wave superconductors with 
very short coherence lengths, our interest is mainly in un- 
derstanding the consequences of magnetic defects when 
there is a clear separation between scales set up by the 
Fermi length kp 1 and the coherence length fo- 

[31] A.V. Balatsky and M.I. Salkola, Phys. Rev. Lett. 76, 
2386 (1996). 

[32] The BCS Hamiltonian transforms as Ho — > — Ho under 
the transformation in which , Fk — ► = ^^k- However, 
because the exchange interaction is invariant under the 
same transformation, it is no longer necessary to have an 
even number of quasiparticle states inside the energy gap 
at nonzero energies. 



14 



